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ABSTRACT 

The observed number density of Lyman-a sources implies a minimum vol- 
ume of the inter-galactic medium that must be ionized, in order to allow the 
Lyman-a photons to escape attenuation. We estimate this volume by assign- 
ing to each Lyman-a emitter the minimum Stromgren sphere that would allow 
half its Lyman-a photons to escape. This implies a lower limit to ionized gas 
volume fraction of 20-50% at redshift z=6.5. This is a lower limit in two ways: 
First, we conservatively assume that the Lyman-a sources seen (at a relatively 
bright flux limit) are the only ones present; and second, we assume the small- 
est Stromgren sphere volume that will allow the photons to escape. This limit 
is completely independent of what ionizing photon sources produced the bub- 
bles. Deeper Lyman-a surveys are possible with present technology, and can 
strengthen these limits by detecting a higher density of Lyman-a galaxies. 

Subject headings: galaxies: high- redshift, (galaxies:) intergalactic medium 



1. Introduction 



The epoch of reionization marks a phase transition in the universe, when the inter- 
galactic medium was ionized. Recent observations oi z > Q quasars show a Gunn-Peterson 
trough (Gunn & Peterson 1965) implying that the reionization of intergalactic hydrogen was 
not complete until z ^ Q (Becker et al 2001, Fan et al 2002). Yet microwave background 
observations imply substantial ionization as early as 2 > 15 (Spergel et al 2003; Kogut et 
al 2003). Transmitted flux seen in some quasar spectra also implies inhomogenities in the 
ionized gas (Oh & Furlanetto 2005; White et al. 2005). These results can be reconciled if 
reionization occurred twice (e.g.. Gen 2003), slowly (e.g. Gnedin 2004) or was substantially 
inhomogenous (Malhotra et al. 2005, Oh & Furlanetto 2005, Gen 2005). Information on the 
state of the inter-galactic gas — the ionized fraction, and spatial distribution of ionized gas — 
is needed, but is scant. 
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Lyman-o; emitters as tests of ionized IGM Lyman-o; emitting galaxies provide an- 
other tool for probing reionization, which is independent of and complementary to of the 
CMBR and Gunn-Peterson diagnostics. Lyman-a visibility tests offer a local probe of neutral 
fractions Xhi ~ 30% (Malhotra & Rhoads 2004 (MR04)), while the Gunn-Peterson trough 
saturates at xhi ~ l%,(Fan et al 2002), and the CMBR polarization provides an integral 
constraint of the ionized gas along the line of sight. 

Because Lyman-a photons are resonantly scattered by neutral hydrogen, Lyman-a line 
fluxes are attenuated for sources in a significantly neutral {{xhi) ~ 0.1) intergalactic medium 
(IGM) (Miralda-Escudc 1998, Loeb & Rybicki 1999, Haiman & Spaans 1999). To zeroth 
order this should produce a decrease in Lyman-a galaxy counts at redshifts beyond the 
end of hydrogen reionization (Rhoads & Malhotra 2001 (RMOl)). This does not mean that 
Lyman-CK emitters will suddenly become invisible at some redshift, but rather that Lyman-a 
flux is strongly attenuated. We consider three physical effects that modify the attenuation 
of the Lyman-a flux. 

Effect of ionized bubbles Each galaxy creates a local bubble of ionized gas. If this is 
large enough, Lyman-a photons are redshifted by the time they reach the neutral boundary 
and thus can escape. Consider a galaxy in a Stromgren sphere of radius /?.,.,, surrounded by 
a fully neutral IGM. The line center optical depth due to scattering by the damping wing 
of neutral gas outside the galaxy's Stromgren sphere is given by r = 1.2pMpc/i?ss (e.g., 
RMOl )} Thus, transmission becomes significant when the physical radius of the HII region 
is > 1.2pMpc (correpsonding to > 9cMpc at redshift z = 6.5), 

For wavelengths separated from line center by AAg^, the optical depth due to neutral 
IGM damping wing is 



where c is lightspeed and where the Hubble constant H ^ 760 km s ^ Mpc ^ x [{1 + z)/7.5]^^'^ 
for a flat cosmology with Hq = 71kms^^ Mpc^\ flm = 0.27, and fl\ = 0.73 (see Spergel et 
al 2003). 

The size of a typical Stromgren sphere that an isolated galaxy creates in the IGM can 
be directly related to its Lyman-a luminosity Liy^ = 10'^^L43 ergs""^, its age t = lO^tgyr, 
and the fraction /esc of its ionizing photons that escape the galaxy to ionize the surrounding 



^In this paper we will write physical Mega-parsecs as pMpc, and comoving Mega-parsecs as cMpc. 
IpMpc = (1 + ^;)cMpc. 




(1) 
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IGM (thereby becoming unavailable for Lyman-a line production) (RMOl). The result, 
conservatively ignoring recombination, is 



Stellar population models for narrowband-selected Lyman-a emitting galaxies require 
ts ^ 1 (and usually ^ 1) to produce the observed range of Lyman-a line equivalent widths 
(Malhotra & Rhoads 2002, MR02). The luminosity function for Lyman-a galaxies at z ~ 6 
shows L^: ~ 10^^'^ erg s"-*^ (MR04), so that galaxies with L43 ^ 1 are rare. Inserting ^ 1 
and L43 < 1 into equation 2 and comparing the result to equation 1, wc sec that r > 2 at 
line center even for reasonably bright Lyman-a galaxies, and that these galaxies' Stromgrcn 
spheres are not sufficient to prevent Lyman-ct line suppression by an order of magnitude or 
more if the surrounding IGM is neutral. 

Velocity Offsets The observed Lyman-a line is usually asymmetric (e.g. Rhoads et al 
2003, Dawson et al 2004) and offset to the red compared to other lines (e.g., Shapley et al. 
2003). This offset is likely caused by absorption or scattering of Lyman-a photons in the blue 
wing of the emission line in the interstellar medium of the emitting galaxy, combined with 
gas motions (or outflows) in that medium. Resonant scattering by the surrounding IGM 
will further suppress the blue wing of the line, and accentuate the asymmetry. The typical 
observed velocity offsets for LBGs with Lyman-a emission is ~ 360km s~^, although the 
velocity offsets decrease for galaxies with higher equivalent widths in the Lyman-a emission 
(Shapley et al. 2003). Such an offset can be incorporated in equation 1 using AAem = ^o^w/c, 
where (the "wind" velocity) denotes the offset of the Lyman-a line relative to the systemic 
velocity of the emitting galaxy. 

If we describe the intrinsic profile as /a(Ao + SX), where Aq = 1215. 67A is the rest 
Lyman-a wavelength, then the unabsorbed line flux is Fq = /^/a(Ao + S\)d{6\). The 
absorbed line flux is F = exp[— r((5A)]/A(Ao -t- 6X)d{6X), and the transmission factor is 
T = F/Fq, which is calculated by integrating over the fine profile. 

Santos (2004) has calculated the transmission factor for a large range of model Lyman-a 
line properties and IGM models His results show that T < 1/3 for a wide range of plausible 
models, even including those with Vyj = 360km s^^ Models with f ^ = invariably show 
T <^ 1/3, unless the IGM is largly ionized. Moreover, the variation of T with Vy^ is fairly 
slow for Vu, > 360km s""^; that is, a large increase of (vu,) would be required to substantially 
raise (T). 
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Galaxy clustering Even though a single Lyman-a galaxy does not produce a large enough 
Stromgren sphere, there is a possibility that fainter galaxies clustered around the observed 
Lyman-a galaxy would contribute enough photons to make the Hll regions around the 
cluster reach the critical size. This would require roughly 10-60 times more photons than 
produced by a typical Lyman-a galaxies observed (Wyithe & Loeb 2005, Haiman & Ccn 
2005, Furlenatto et al. 2004) So far the evidence for such clustering is mixed: Malhotra et 
al. 2005 see clustering in the Hubble Ultra Deep field which was chosen to include a 25th 
magnitude galaxy at z=5.8, so do StiaveUi et al 2005 around a SDSS quasar at z=6.2; but 
Rhoads et al. (in prep) see no significant excess around an actual z=6.5 Lyman-a galaxy. 

We here introduce a new version of the Lyman-a rcionization test wherein aach observed 
Lyman-a galaxy is seen because the IGM near that galaxy is ionized in an otherwise neutral 
medium. Each Lyman-a galaxy thus implies the presence of a certain volume of ionized IGM. 
We combine this ionized volume per source with the observed number density of Lyman-a 
galaxies to obtain a lower bound on the volume ionized fraction of the IGM. This limit is 
completely independent of whether the ionizing photons come from the Lyman-a galaxy, its 
unseen neighbors, or any other source. 



2. Ionized volume estimates 

Calculating the Transmission We need to calculate the Lyman-a flux transmission 
expected from a galaxy with a given Stromgren sphere radius, Lyman-a line velocity offset, 
and IGM neutral fraction. To do so, we must make some reasonable assumptions about the 
line width and shape. Observed Lyman-a emission lines from high redshift galaxies are well 
described by a truncated Gaussian profile (Hu et al 2004, Rhoads et al 2004): 

f .^E^exnf-^) iiSXX) 
fx{Xo + SX) = { 2.^; ''^^-"^ (3) 

( if < 

where Aq = 1215.67(1 -|- v^/c)A is the Lyman-a central wavelength in the frame of the 
IGM surrounding the galaxy, and the unabsorbcd line flux is Fq. The absorbed line flux is 
F = cxp[— r(5A)]/A(Ao + 6X)d{6X), and the transmission factor is T = F/Fq. We set 
a = 1.49A, or rest frame FWHM of 1.75A (from the truncation point Aq to the half-peak 
point on the red side of the line) - chosen to match typical observed high redshift Lyman-a 
fines. 

We determined T by numerical integration for a grid of and Rss- We then inverted 
the result to determine the Stromgren sphere radii corresponding to 25%, 33%, 50%, and 
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70% transmission for a range of Lyman-a velocity offsets Vy,- The result is summarized in 
figure 1. 

The Minimum Ionized Fraction: Analytic results We next combine our calculated 
minimum Stromgren sphere radii with the observed number density of z ~ 6.5 Lyman-a 
emitting galaxies to place a lower bound on the volume ionized fraction of the intergalactic 
medium. The problem is closely related to the void probability function, or to a counts-in- 
cells formalism. The volume neutral fraction xv,n = 1 — xv,ion can be found by calculating 
the probability that an arbitrary point in space lies in a void of radius > Rgs in the Lyman-a 
galaxy distribution, or (equivalently) by calculating the probability that a randomly placed, 
spherical cell of radius Rgs will contain no Lyman-a galaxies. 

The zero order estimate of the ionized volume fraction is simply xy^ion ^ T = nVss- Here 
n is the volume number density of the detected Lyman-a emitting galaxies, = (47r/3)i?2g 
is the volume of the minimum size Stromgren sphere that renders them visible, and J- is the 
fiUing factor of these spheres. The estimate xv,ion ~ ^ is accurate for xv,ion <S 1. At larger 
ionized fractions, we need to correct for overlap of the Stromgren spheres. 

If the locations of the Lyman-a galaxies are uncorrelated, the residual neutral volume 
fraction becomes xy^n = exp(— JF) and the corresponding ionized fraction is xy^ion = 1 — 



The degree of overlap will be further enhanced by spatial correlations amonff galax- 
ies. To account for this overlap, we can assume that the observed Lyman-a emitters are 
drawn from a distribution with two-point angular correlation function ^{r) ^ (^/^o)^- 
We take values ro ~ 4cMpc and 7 = 1.8, which are typical of high redshift galaxy popu- 
lations and in particular are consistent with observed Lyman-a galaxy populations (e.g., 
Ouchi et al. 2003, Kovac et al. in prep.). This correlation function implies that the 
mean number of neighbors within distance r of an arbitrarily selected galaxy should be 
{Np) = (47r/3)r3n (1 + 3(r/ro)-^/(3 - 7)) (Peebles 1980, eq. 31.8). Given a pair of ionized 
spheres with radius Rss, whose centers are separated by distance r (with < r/Rgg < 2), 
the volume of their overlap region is Vss x (1 — 0.75{r / Rss) + {r / Rs s)^ / 16) ■ We integrate the 
expected volume ''lost"' to overlap within a survey by combining the two foregoing results: 



While this calculation accounts for correlations, it omits third order and higher terms in J^. 
The coefficient for the J-'^ term would depend on the three point correlation function. Thus 
equation 4 will tend to underestimate the volume ionized fraction, and is conservative. 



exp(-.^) = Er=i(-l)^-'-^Vj!- 
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Fig. 1. — This figure shows the minimum radius of the ionized bubble needed for 70% (dot 
dashed hne), 50% (sohd hne), 33% (short dashed hne) or 25% (long dashed line) of Lyman-a 
line flux to escape as a function of the velocity offset between the Lyman-a emission line 
and the intrinsic velocity of the galaxy. Currently, a transmission factor of less than 50% 
is ruled out by MR04, and the best guess for velocity offset is ^ 350km s~^; which would 
imply Rss — 0.93 physical or 7 comoving Mpc. 
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Minimum Ionized Fraction: Numerical simulations To achieve accurate results at 
higher source densities, we simulated correlated distributions of galaxies and directly calcu- 
lated the volume fraction enclosed by the union of their Stromgren spheres. We performed 
these calculations for a range of correlation lengths, Stromgren sphere sizes, and source den- 
sitites. We generated the set of correlated {x, y, z) triples using a Mandelbrot-Levy random 
walk prescription (see Peebles 1980, section 62). A power law distribution of step lengths 
s, p{s) oc s^^'^, ensures a power law correlation function with the desired slope 7 = —1.8. 
The correlation length tq is related to the minimum step length (which is always <S Rgg and 
<^ To) and to the number of independent random chains occupying the sample volume. To 
tune the correlation length to a desired value, we vary the number of chains used and also 
(for small correlation lengths) add a suitable number of entirely uncorrelated (x, z) triples. 

We verified that this prescription generated the desired two-point correlation properties. 
The volume used for the simulations was ~ lO^cMpc^, and the volume fraction was calculated 
on a grid with cell size = -Rss/2. (A smaller cell size does not significantly change the results) 
Finally, wc verified that the simulations reproduce the analytic results for an uncorrelated 
distribution at any filling factor, and for a correlated distribution at < 1 (where the 
and higher terms are small). 

Results of the simulations are shown in figure 2, where we show curves of 
function of the Lyman-a galaxy number density for a plausible range of Rgg and tq values. 
Were we to plot xy^ion as a function of JF (rather than n), wc would sec that Xv^ion is a function 
of just two variables, T and the ratio {vq/Rss), as one would expect from the discussion in 
section 2. 

Effect of a Partly Ionized Ambient Medium Up to now, we have considered Strom- 
gren spheres embedded in a fully neutral ambient intergalactic medium. If this medium is 
already partially ionized outside the bubbles, then Lyman-ct fiux is transmitted more eas- 
ily, and the Stromgren spheres required to achieve a particular Lyman-a transmission are 
smaller. However, this scenario requires a large amount of ionized gas mixed into the ambi- 
ent medium. Figure 3 shows the net ionized fraction of ionized gas (counting the Stromgren 
spheres -|- ambient ionized gas) as a function of ambient neutral fraction for — and 
360km s-^ 
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Fig. 2. — The volume ionized fraction as a function of observed Lyman-a source density 
is plotted for a range of minimum Stromgren sphere sizes and correlation lengths. From 
highest to lowest, these sets of curves correspond to Rss = S.OpMpc (yellow), 2.5pMpc 
(red), 2.0pMpc (magenta), 1.4pMpc (green), 1.2pMpc (cyan), l.OpMpc (black) and O.SpMpc 
(blue). For each set of curves, the hne style indicates the correlation length: tq = (dotted), 
and 4cMpc (solid). Vertical lines at 1.2 x 10"'* and 1.8 x 10~^cMpc~^ mark the range of 
Lyman-a number density observed by Taniguchi et al (2004). For the best-guess parameters 
1 < Rss < 1.4pMpc and tq ^ 4cMpc, the Taniguchi data imply a lower bound of 0.2 to 0.5 

on Xv,ion- 
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Ambient neutral fraction 

Fig. 3. — The ionized gas fraction in the IGM is plotted as a function of the neutral fraction 
in the ambient medium (outside Stromgren spheres). We have assumed a number density 
of 1.5 X 10~^cMpc~^ (Taniguchi et al 2004). Two solid (black) curves show the result for 
Vw = 0. The lower curve shows the amount of ionized gas enclosed in Stromgren spheres 
around Lyman-a galaxies, and the upper curve shows the total ionized gas fraction (including 
that in the ambient medium). The dashed (blue) curves show the same quantities for = 
360km s~^. 
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3. Summciry and Future outlook 

Observations of Lyman-a; emitters provide a powerful method to probe the neutral 
fraction of the IGM and the volume of the ionized gas. In an earlier paper we explored the 
effect on the Lyman-a luminosity function and constrained the average neutral fraction of 
the IGM to be (a;„) < 30-50% (MR04, see also Stern et al. 2005). Further investigations by 
Haiman & Oh 2005 and Furlennato et al. 2005 agreed with these estimates. In this letter we 
place a lower limit on the volume in ionized IGM xv,ion — 20 — 50% for best known parameters 
of transmission factor T — 0.5, source density n{lya) — 10~^, and v-w — 360 — 0km s~^. 

With better observations, these limits can be substantially improved. Gurrently, the 
largest Lyman-a galaxy sample at z—6.5 comes from Subaru observations (Taniguchi et al. 
2005) which reach just about L*. Going deeper is then expected to yield a higher density 
of Lyman-a emitters. If the current observed number density of sources could be increased 
by a factor of 2 or 3 with deeper observations, the lower limit to Xv^ion would increase by a 
factor of 2. Bigger samples (both at z=6.5 and z=5.7) would also improve the estimate of 
the attenuation factor. MR04 concluded that attenuation of a factor of two was marginally 
allowed by present data. If we could reduce this limit to 1.4, i.e a transmission factor of 70%, 
we would conclude that the minimum Stromgren sphere radius needed would be 3 pMpc for 
offset velocity of zero, and 2.6 pMpc for a velocity offset of 350 km s~^ (Figure 1). In that 
case the Taniguchi et al. 2005 data would imply nearly xv,ion — 100%. Moreover, we should 
not forget that we are using the smallest value of Rss that allows the transmission of the 
Lyman-a line. The ionized bubbles could certainly be larger. 

Measurement of the velocity offset between Lyman-a redshift and systematic velocity 
of the galaxy will allow us to further refine the Stromgren sphere radius (see Figure 1). 
Gurrently our estimates come from such measurements on Lyman break galaxies (LBGs), 
for which the higher the EW of the Lyman-a line, the smaller the velocity offset (Shapley et 
al. 2003). But the overlap between the LBG samples and the Lyman-a samples is minimal; 
only the LBGs with highest EW (EW ~ 50 A) would even be selected by the Lyman-a 
searches using Narrow-bands. So it is plausible that the velocity offsets for the Lyman-a 
emitters (median EW ^ lOOA; (MR02, Dawson et al. 2004) would be even lower. 

Spectroscopy of large, complete samples would also yield the information on their dis- 
tribution and clustering in the line-of-sight dimension, which can be used to determine the 
overlap in the ionized bubble around the galaxies. Then we can go beyond the zeroth order 
estimate of xv,ion to learn about the topology of the ionized bubbles and thus determine 
whether substantial overlap has happened (Rhoads, in prep). Pre-overlap phases of reioniza- 
tion can also be detected by heightened correlation between Lyman-a sources (Furlanetto et 
al. 2005). If the sources required for ionizing photons are strongly clustered on f« lOOcMpc 
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scales as observed (Malhotra et al. 2005) or predicted (Cen 2005), we should see dramatic 
field-to-field variations in Lyman-o; number density prior to reionization. 

These tests only require the detection of statistically useful samples of Lyman-a emit- 
ters and are eminently practical with present technology at optical wavelengths, and their 
extension to redshifts z > 1 through near-infrared observations is close at hand. 

We thank Michael Strauss for useful discussions. 
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